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ABSTRACT 

We present direct evidence for anomalous microwave emission in the Perseus molecular cloud, which 
shows a clear rising spectrum from 11 to 17 GHz in the data of the COSMOSOMAS experiment. By 
extending the frequency coverage using WMAP maps convolved with the COSMOSOMAS scanning 
pattern we reveal a peak flux density of 42 ± 4 Jy at 22 GHz integrated over an extended area of 
1.65° x 1.0° centered on RA=55.4° ± 0.1° and Dec=31.8° ± 0.1° (J2000). The flux density that we 
measure at this frequency is nearly an order of magnitude higher than can be explained in terms of 
normal galactic emission processes (synchrotron, free-free and thermal dust). An extended IRAS dust 
feature G159.6-18.5 is found near this position and no bright unresolved source which could be an 
ultracompact H n region or gigahertz peaked source could be found. An adequate fit for the spectral 
density distribution can be achieved from 10 to 50 GHz by including a very significant contribution 
from electric dipole emission from small spinning dust grains. 

Subject headings: diffuse radiation dust, extinction — ISM; individual (G159.6-18.5) — radiation 
mechanisms: general — radio continuum: ISM 


1. INTRODUCTION 

The search for ‘anomalous microwave emission’ or 
‘dust-correlated microwave emission’ was started by the 
surprising statistical correlation of COBE DMR obser¬ 
vations at centimetri c wavelengths with DIRB E dust 
emissi on at 140^m llKogut et all 1 19961) . iLeitch et all 
(I1997TI found a similar correlation with the OVRO 
ring observations at 14.5 and 32 GHz. Ordinary 
thermal emission from dust at these relatively low 
frequencies is expected to be orders of magnitude lower 
than the intrinsic CMB fluctuations. It was initially 
suggested that this anomalous emission had its origin 
with free-free processes associated with the warm 
ionized interstellar medium. The lack of the expected 
Ha from this ionized gas would r equire tempe ratures 
hotter than 10 6 K. However IDra.ine and La.zarian (Il998f) 
argued against such high temperature on energetic 
grounds and proposed an alternative mechanism 
based on electric dipole emission from small dust 
grains. Further evidence for anomalous microwave 
emission has been found in data from other CMB 
exper iments such as Saskatoo n et, a.l .1 

119971). 19GHz iT le Olivei ra-Costa et^TTTTp D . Tenerife 
(Ide Oliv eira-Costa et, alJ I1999L 120021 120 1)411. Python 
V llMuk herj^ e t, alJ l2003lb WMAP llLagachel 2003; 
Fiukh einerl 120041) and the Green Bank Galactic Plane 
Survey (GBGPS) (IFinkbeiner. Langston and Minted 
120048b 

Despite this growing evidence, the exact nature and 
distribution of this new foreground remains elusive 
and part of this is because there is no compara¬ 
ble survey below 20 GHz to compl e ment the WMAP 
data. For instance, iBennett et al.1 : 2003 a: find that 
WMAP foregrounds can be fitt ed by a pro posed flat¬ 
ter component of synchrotron. IFinkbeiner! 0 2004 1 and 


Ide Oliveira-Costa et all 120041) point out that Tener¬ 
ife and GBGPS data provides evidence for rising 
spectrum at 10 and 15 GHz which is incompatible 
with synchrotron e mission. A similar s urvey to the 
Tenerife experiment iGutierrez et a DEqqcJ in frequency 
and sky coverage but with degree scale resolution is 
now available with the COSMOSOMAS experiment 
(IFernandez-Cerezo et, all 120051 iHildebrandt, et all 120051) 
of the Instituto de Astrofisica de Canarias (IAC). 

The COSMOSOMAS experiment consists of two cir¬ 
cular scanning instruments operating at the Teide ob¬ 
servatory (altitude 2400 m, Tenerife). They have pro¬ 
duced 0.8-1.1° resolution maps of ~ 10000 square 
degrees in the sky in four frequency bands (centered 
at 10.9, 12.7, 14.7, 16.3 GHz). The average sensi¬ 
tivity of each map is in the range 80 — 120/zK per 
beam. The b asic o bserving strategy is described in 
iGallegos et al ] lt200ilh Details of the survey by the two 
instruments where a signifi cant du st co r relation i s found 
at h > 20° is giv en in I Fernandez-Cerezo et ill. Hill); 
Hilde brandt et all lj2005l) . In this letter we present re¬ 
sults of a search for sites of dust-correlated emission with 
a rising flux density spectrum in the first release of the 
COSMOSOMAS maps. 

2. THE SEARCH FOR AREAS OF RISING SPECTRUM 
WITH COSMOSOMAS 

Since we are focusing on bright discrete sources we 
cho se to generate a list of possible sources with sextrac- 
tor iBertin and Arnoutslll996 f using a detection thresh¬ 
old of 4 sigma on the more sensitive 13 GHz COSMOSO¬ 
MAS map (+22.0° <Dec< +45.0°). Since in the Auriga- 
Perseus galactic crossing the confusion is lower and the 
sources are more clearly separated, we have centered our 
search for sources of rising spectrum in this region which 
is shown in fig 1. 
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Fig. 1.— Comparison of various survey maps over the Auriga-Perseus region. 


A few of the objects turned out to be artifacts on close 
inspection and fluxes were computed for all remaining 
true sources in each of the COSMOSOMAS frequency 
bands using a 3 by 3 pixel box (1°RA by l°Dec). One 
source (RA = 55.4 ± 0.1°, Dec = +31.8 ± 0.1° J2000) 
stood out with a clear steep rising spectrum between 11 
and 17 GHz with an index of a = +1.4 ± 0.3 (S„ oc v a ). 
This source lies in the Perseus OB2 area at an interme¬ 
diate latitude b=-18° and is also seen as a bright source 
in the 100/im DIRBE map which appears to correspond 
to the extended IRAS source G159.6-18.5. It also ap¬ 
pears to be extended in the COSMOSOMAS maps with 
a FWHM w 2.0°. The average temperatures of this pixel 
box are shown in Table 1 showing the almost constant 
thermal spectrum over the COSMOSOMAS fr equencies. 

To m ake a comparison, the WMAP data lBeunett et afl 
(l2003al) have been degraded to the common COSMOSO¬ 
MAS resolution (set at 1.12° ), convolved with the simu¬ 
lated circular lock-in analysis and pixelated in a regular 
1/3° Cartesian grid in RA and Dec. The lockin analysis 
includes the removal of the first 7 harmonics to suppress 
the 1/f noise from receiver and atmosphere. The effect 
on the WMAP maps is shown by comparing the original 
map fig.1(a) and the convolved map fig.1(b), where the 
loss of broad angular features is evident with the creation 
of negative lobes north and south of bright objects in the 
center of the map. This figure also gives a visual exam¬ 
ple of this dust correlation through the strikingly similar 
WMAP K (22 GHz) (fig.la) and DIRBE lOOpm(fig.lc) 
maps in the galactic Auriga-Perseus region. 

Two of the detected sources in this area are bright 
diffuse Hu regions, NGC 1499 and IC4 10 that we use as 
a cal ibration check. Using 1420 MHz llReich and R.eichl 
If 988 1 and the five WMAP template maps we verified that 


the emission from these objects has a spectral index of 
—0.15 ± .02 and —0.20 ± .03, respectively. These indices 
are consistent with that expected for free-free emission 
of -0.12 between 11 and 17 GHz. For both sources we 
find the measured COSMOSOMAS fluxes are generally 
lower by « 6% than the predicted values from the fits. 
This is consistent with the absolute calibration error of 
the COSMOSOMAS maps with is better than 10% in 
the Auriga-Perseus region. 

In order to get the dust correlated emissivity (fj ,K 
[MJy] -1 sr) for G159.6-18.5 and to quantify the anoma¬ 
lous microwave emission we have carried out a correlation 
with the 100/mr DIRBE map fig. 1(c) passed through the 
COSMOSOMAS scanning. We chose a 11 by 11 pixel 
box, an area of 3.1° x 3.7°, centered on the pixel nearest 
to the center of this source. This covers the bulk of the 
emission allowing a good range of pixel values without in¬ 
cluding much background. The correlation is calculated 
by the linear fit of temperature (/iK) of the COSMO¬ 
SOMAS and WMAP pixels to the intensity (MJy/sr) of 
the DIRBE pixels using a standard least squares method. 
Table 1 gives the values for the slopes and errors showing 
a tight correlation with DIRBE. The equivalent analysis 
with the 408 and 1420 MHz templates suggests free-free 
emission is present but with less significance since this 
emission appears more compact and displaced with re¬ 
spect to the DIRBE emission. 

3. G159.6-18.5 AND ITS SPECTRAL ENERGY 
DISTRIBUTION 

G159.6-18.5 lies within the Perseus Molec¬ 
ular clou d complex at a di stanc e of 260 pc 
llCernicharo. Bachiller and Duvertl[l98 5). It appears in 
IRAS maps as a slightly broken ring with a diameter 
of « 1.5° composed of several diffuse knots around the 
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TABLE 1 

Spectrum of G159.6-18.5 


V 

(GHz) 

Temperature 

(mK) 

Correlation 
(/rK MJy— 1 sr) 

Flux 

(Jy) 

0.408 

1346 ± 214 

56000 ± 70000 

6.3 ± 7.8 

1.420 

144 ± 14 

4260 ± 1410 

7.3 ± 2.0 

10.9 

2.264 ± 0.096 

59.6 ± 0.8 

17.0 ± 0.1 

12.7 

1.798 ± 0.098 

45.1 ± 1.9 

18.7 ± 0.4 

14.7 

1.839 ± 0.072 

44.5 ± 2.0 

26.2 ± 0.6 

16.3 

1.682 ± 0.065 

41.9 ± 4.0 

32.6 ± 1.5 

22.8 

1.195 ± 0.035 

33.9 ± 0.3 

42.3 ± 0.2 

33.0 

0.547 ± 0.016 

15.7 ± 0.3 

40.3 ± 0.4 

40.9 

0.296 ± 0.009 

7.7 ± 0.3 

33.9 ± 0.7 

61.3 

0.134 ± 0.004 

3.3 ± 0.3 

34.7 ± 1.8 

93.8 

0.116 ± 0.003 

2.7 ± 0.4 

77.5 ± 4.3 

1250. 



96800 ± 700. 

2143. 


131000 ± 1250. 

3000. 



64400 ± 100. 


rim and a couple in the center, all with intensities of 
100-200 MJy/sr at 100/irn. These knots are still bright 
at 12/irn with intensities of 5-10 MJy/sr suggesting the 
presence of transiently heated small dust grains. 

At first it was thought to be a supernova remnant 
(SNR) bv lFiedler et al]1119941) . but it was later argued 
to be a H II region driv en by the 09.5-BO V star HD 
278942 (Ide Zeeuw et al.lll999li found to be at the center 
of the ring. lAndersson et a.U nu have made an ex¬ 
tensive overview of the nature of this dust ring including 
interferometric observations at 408 and 1420 MHz and 
conclude the ring is a ruptured blister HII region which 
has emerged from the outer edge of the cloud. Their 
flux density estimates are F(A08MHz) = 2.0+q ® Jy and 
F(1420MHz) = 1.41q’^ Jy allowing for incomplete sam¬ 
pling of the u — v plane. This gives a spectral index of 
a = —0.3/li'i, which is consistent with an optically thin 
Hu region (a = —0.1). They also found many radio 
sources, but all are steep spectrum apart from the flat 
spectrum quasar 4C+32.14, which appears at the NW 
edge of the ring. This was removed from all the maps 
under consideration here assuming a flux of 2 Jy at all 
frequencies. 

To investigate the emission mechanisms occurring in 
this object we need to estimate its spectral energy dis¬ 
tribution (SED) over a wide range of frequencies. Since 
flux on large angular scales is lost due to the COSMO- 
SOMAS lock-in analysis we are forced to fit the data to 
a convolution of a sky density flux model. We adopted a 
simple elliptical Gaussian for G159.6-18.5 based on a fit 
to the original WMAP Ka map. We obtained a FWHM 
ellipse for the unconvolved model of 1.6° x 1.0° PA = 
51°. The expected beam spread function for COSMO- 
SOMAS was then fitted to all the template maps using a 
least squares method allowing for a zero level offset. For 
these fits we obtained the amplitudes and subsequently 
the fluxes, which are listed in table 1 and plotted in fig.2. 
The two low frequency surveys do not show the extended 
emission, but unresolved emission from the center of the 
ring so we used the sky model corresponding to a point 
source. The error in the estimate of the FWHM of the 
sky model of 0.1° leads to an overall absolute calibration 
error of 10% although the relative calibration between 
frequencies is better than 2%. 


4. LIMITS ON STANDARD EMISSION MECHANISMS 

The main galactic microwave emission mechanisms are 
synchrotron, free-free and thermal dust emission. Syn¬ 
chrotron can be estimated from the low frequency sur¬ 
veys, free-free from Ha emission and thermal dust from 
DIRBE. 

The values from the 408 and 1420 MHz surveys for 
G159.6-18.5 suggest no significant synchrotron emission 
and some free-free emission at the these frequencies, 
consistent with the flu x density measurements made by 
lAndersson et al.l ( 2000 1 inside the ring. A “bean” shaped 
patch of Ha emission inside the dust ring is detected in 
this region by t he Wisconsin H Al pha Mapper project 
(WHAM) ( lHaffner. Reynolds and Tuftell999l) . This sug¬ 
gests the Ha emission is due to diffuse Hii and is shin¬ 
ing through a thinner or ruptured part of the dust 
shell. Estimating free-free emission from the Ha emis¬ 
sion is difficult due to the high extinction which also 
varys strongly across the ring and the unknown frac¬ 
tion of dust in front of the H II region, but we only re¬ 
quire a rough indication o f the level of em i ssion. We 
adapted the method that IDickinson et alJ ( 2003 1 de¬ 
scribe to estimate galactic free-free emission from Ha. 
Assuming an extinction factor Rv = 4.0 as expected 
for molecular clouds we use an expression for Ha ex¬ 
tinction A( Ha) = 3.33 E (B- V) together wit h the E(B- 
V) maps of iSchlegel. Finkbeiner and Davisi 1998) and 
find a mean corrected Ha intensity over the ring of 120 
Rayleighs. At 22 GHz, close to the peak source emission, 
the free-free brightness temperature to Ha intensity ra¬ 
tio (T^/Iho), assuming electron temperature 8000 K, 
is 12.3/iAT/Rayleighs. Integrating over the area of the 
ring this corresponds to a flux density of 4.5 Jy, similar 
to the 5.6 Jy estimated by extrapolating from the 1420 
MHz flux density derived from the Reich and Reich map. 
This is nearly an order of magnitude below the estimated 
flux density of 42 Jy from the WMAP K data. 

Thermal dust emission can be fitted to the DIRBE and 
WMAP W points as shown in fig.2, with a dust temper¬ 
ature of 19 K with opacity at 100/.im(noo) of 6.0 x 10 -3 
and emissivity index 1.55. This is what is expected for 
dust in a “warm neutral medium” (WNM). This WNM 
region might be generated by interaction with the nearby 
B0 star. A v ery co ld dust grain model as tried by 
iCasassus et a.U ( 2004 ) with CBI observations of anoma¬ 
lous microwave emission from the Helix nebula would not 
be consistent with COSMOSOMAS and WMAP points. 

There are two known types of object which show a 
convex spectrum at microwave frequencies; ultracompact 
Hii regions (UCHii) and Gigahertz Peaked Spectrum 
(GPS) sources which we assess as possible sources to pro¬ 
vide the excess flux density at 22 GHz. 

UCH II regions produce free-free emission which is self- 
absorbed at low frequencies due to extremely high densi¬ 
ties and the high d ust opacity can hide the Ha emission. 
Indeed iMcCullough and Ghenl 12002 1 suggested that an 
UCH II could exp l ain the microwave emission detected by 
iFinkbeiner et a ,1. (2002) in the diffuse Hu region LPH 
201.663+1.643, but it did not seem to have the neces¬ 
sary FIR luminosity and so anomalous emission was not 
ruled out. In the case of the Perseus cloud its distance is 
known and UCH n regions being amongst the brightest 
objects at 100/im would have a flux density in excess of 
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Fig. 2.— Spectrum of G159.6-18.5. Points are shown for COSMOSOMAS, WMAP , DIRBE and low frequency (408 & 1420 MHz) surveys 
convolved to give COSMOSOMAS equivalent points. The size of the symbols correspond to the sky model uncertainty of 10%. Estimates 
are made of three foregrounds; free-free, rotational dust and vibrational dust. Vibration dust corresponds to T^ust = 19.OK and emissivity 
index 1.55. A Toy dust model similar to the one used by Finkbeiner is used which consists of a linear combination of Draine and Lazarian 
models for Warm Neutral Medium and Molecular Cloud (0.8*WNM + 0.3*MC). 


1 MJy. None of the IRAS point sources within l°of the 
COS MO S OMA S source center hav e IRAS colors which 
meet IWoorl and Churchwell ( 198 9) criterion for UCHll 
regions and are all less than a few hundred Janskys. Such 
a source could be behind the cloud at a distance of > 10 
kpc, but at galactic latitude -18°this seems very unlikely. 

GPS sources are high red-shift radio sources in which 
the radio jets have been highly confined and the syn¬ 
chrotron emission is self-absorbed. The majority of GPS 
sources have a peak around 1 GHz, b ut there are a 
few rare ones such as 1459.9+3337 (lEdge et all 1 1 9981) 
which peak at 30 GHz. A GPS source having an in¬ 
dex of +1.4 on the low frequency side as suggested by 
COSMOSOMAS data would have a flux density at 4.85 
GHz of ft! 5 Jy. A searc h in the Green Bank survey 
(GB6) llGregorv et alJll99 6) within 1.5°of the COSMO¬ 
SOMAS source center finds the brightest source, apart 
from 4C+34.14, J0340+3209 which is 0.47°away with a 
flux d ensity of 482 m Jy. A cross check in the NVSS cata- 
logue lCondon et aT . ( Tl998f l over the same area finds this 
source again as the brightest with a flux density of 1.6 
Jy showing it to be a steep source with no significant 
emission at the COSMOSOMAS frequencies. 

The extended nature of the COSMOSOMAS source 
argues against a radio point source and hence against 
either of these objects. 

5. DISCUSSION 

iDraine and Lazarian (Il998t 119991 ) have made detailed 
calculations for the expected spectrum from electric 
dipole emission of very small spinning dust grains and 
magnetic dipole radiation from hot ferromagnetic grains. 


It is still uncertain the exact mass and dipole moment 
distribution these grains will have. 

As a test of the spinning dust hypothesis we can check 
if it can produce the measured spectrum with respect to 
the thermal dust emission. Draine and Lazarian’s pre¬ 
dictions are made in terms of emission in intensity per 
hydrogen column density in units of Jy/sr cm 2 /H. We es¬ 
timate a column number density of hydrogen atoms N(H) 
= 1.3 x 10 22 using th e canonical factor of 2.13 x 10 24 H 
cm~ 2 = unit r 10 o as in iFinkbeiner. Langston and Minted 
lf200l . There is not enough information to make a de¬ 
tailed model of emission taking into account the varying 
conditions throughout the cloud, so we use an integrated 
sum of the WNM model A of Drain and Lazarian for the 
dense region facing the B0 star with a contribution of the 
“molecular cloud” (MC) model for dust warmed further 
in the cloud. Assuming the solid angle of the emission 
is that of the sky model used in section 2 (1.6° x 1.0°) 
we find a qualitatively good agreement with the observa¬ 
tions using 0.8 WNM and 0.3 MC of the predicted curves 
(see figure 2). Although we have not attempted it, there 
maybe magnetic dipole emission models which can fit the 
data also. 

This anomalous signal in this cloud is so evident be¬ 
cause the competing free-free emission is low. If the free- 
free emission was five times brighter, for instance if the 
star was more massive and more embedded, then the 
detection would be marginal. Yet the stars interaction 
seems important to excite the thicker parts of the dust 
shell still embedded in the cloud. It seems unlikely that 
objects such as G159.6-18.5 could be at high galactic lati- 
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tudes, but it may prove to be a relatively close example of 
anomalous emission which can be studied in detail. The 
bright FIR luminosity and associated microwave contin¬ 
uum emission may mean other such objects exist at low 
galactic latitudes but have been classified as UCHii re¬ 
gions. 

6. CONCLUSION 

Our search for rising spectrum sources in the COS- 
MOSOMAS data finds anomalous microwave emission 
associated with the IRAS dust feature G159.6-18.5 in 
the Perseus molecular complex. We measure a rising 
flux spectrum from 11 to 17 GHz with a spectral index 
of +1.4. Extended frequency coverage using cross cal¬ 
ibrated WMAP data reveals the expected spectral roll 
off at both high and low frequencies required for a solid 
detection of anomalous emission. The bulk of the emis¬ 


sion seems to come from the thicker dusty wall on the 
molecular cloud side of the ring rather than the center of 
the ring itself. 

We discard an ultracompact Hu region or a Giga¬ 
hertz Peaked Spectrum source as the cause of this 
microwave emission since they would appear as point 
sources whereas both COSMOSOMAS and WMAP find 
an extended diffuse feature. A simple combination spin¬ 
ning dust models can explain the excess microwave emis¬ 
sion. 

We acknowledge the support of the Instituto de As- 
trofisca de Canarias and the staff of Teide Observatory 
for the construction and operation of the COSMOSO¬ 
MAS experiment. Partial funding was provided by grant 
AYA2001-1657 of the Spanish Ministry of Science and 
Education. 
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